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Abstract. In order to extend the Landauer formulation of quantum transport to correlated fermions, we 
consider a spinless system in which charge carriers interact, connected to two reservoirs by non-interacting 
one-dimensional leads. We show that the mapping of the embedded many-body scatterer onto an effective 
one-body scatterer with interaction-dependent parameters requires to include parts of the attached leads 
where the interacting region induces power law correlations. Physically, this gives a dependence of the 
conductance of a mesoscopic scatterer upon the nature of the used leads which is due to electron interactions 
inside the scatterer. To show this, we consider two identical correlated systems connected by a non- 
interacting lead of length Lc- We demonstrate that the effective one-body transmission of the ensemble 
deviates by an amount A/Lc from the behavior obtained assuming an effective one-body description for 
each element and the combination law of scatterers in series. A is maximum for the interaction strength U 
around which the Luttinger liquid becomes a Mott insulator in the used model, and vanishes when U — » 
and U oo. Analogies with the Kondo problem are pointed out. 

PACS. 71.27.+a Strongly correlated electron systems; heavy fermions - 72.10.-d Theory of electronic 
transport; scattering mechanisms - 73.23.-b Electronic transport in mesoscopic systems 



1 Introduction 

. In Landauer's formulation of quantum transport [Q, the 
measure of the conductance g of a coherent system is 
formulated as a scattering problem between incoherent 
electron reservoirs. In a two-probe geometry, the system 
is connected to two reservoirs via leads. For large elec- 
tron densities, the Coulomb interaction is screened and 
the Coulomb to kinetic energy ratio r s is small. One has 
essentially a non-interacting system of Fermi energy Ep, 
where the occupation of the one-body levels is given by 
a Fermi-Dirac distribution at a temperature T. The sys- 
tem acts as a one-body scatterer and its residual conduc- 
tance g(T — ► 0) is given (in units of 2e 2 /h for single chan- 
nel leads and spin degeneracy) by the probability \t(Ep)\ 2 
of an electron of energy Ep to be elastically transmitted 
through it. 

The problem of describing coherent electronic trans- 
port becomes more complex in the case where the car- 
rier density is low inside the scatterer, the screening ceas- 
ing to be effective and the electrons becoming correlated. 
Such situations occur in quantum point contacts of trans- 
verse size smaller than the Fermi wavelength, where a 
0.7 (2e 2 //i) structure is observed PJ, and can be expected 



in molecules atomic chains or contacts 0], quantum 
dots where few electrons might form a correlated solid, a 
charge density wave, a Mott insulator, etc. In these cases, 
the electrons are transmitted from one Fermi reservoir to 
another through a many-body scatterer. To extend Lan- 
dauer's approach to such systems, at least for low tem- 
peratures and bias voltages, one needs to reduce the bare 
many-body scatterer to an effective one-body scatterer 
with interaction-dependent parameters. This task will be 
hopeless for an isolated system where electrons interact 
with a large interaction strength U, but becomes possible 
when leads where electrons do not interact are attached 
to it. This has been numerically demonstrated in previ- 
ous works jSE] using the embedding method, which al- 
lows to extract 5.6 7 8 9 ilOUl 1| the effective coefficient 
\t(Ep,U)\ 2 from the persistent current of a large non- 
interacting ring embedding the many-body scatterer. Us- 
ing the same method, we show that it is not the region 
where the electrons interact which acts as an effective 
one-body scatterer with renormalized parameters, but a 
larger region where the many-body scatterer induces cor- 
relations. This problem is somewhat similar to the Kondo 
problem, which can be solved using Wilson's numerical 
renormalization group (NRG) Instead of using the 
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Fig. 1. Scheme of the ring pierced by a flux (j> used for the 
embedding method. The correlation cloud induced by Ls in- 
teracting sites upon the auxiliary lead is sketched in grey. 




Fig. 2. Scheme of the set-up with two identical many-body 
scatterers connected by Lc sites where the carriers do not in- 
teract. When Lc is small, the two correlation clouds sketched 
in grey overlap and the effective one-body scatterer of the en- 
semble is not given by the effective one-body scatterer of each 
element and the combination law of one-body scatterers in se- 
ries. 



NRG method, we use the density matrix renormalization 
group (DMRG) method 13, l4j for a non-interacting ring 
embedding the many-body scatterer, as sketched in Fig.^ 
The modulus of the effective one-body transmission am- 
plitude \t(Ep, U)\ is obtained from the persistent current 
of the ring extrapolated to infinite lead length, while its 
phase a is given by the Friedel sum rule. 

After a study of the contained extrapolation, we ap- 
ply the embedding method to determine the effective 
total transmission coefficient ^(Lp, U)\ 2 of two identi- 
cal many-body scatterers in series, connected by a non- 
interacting lead of size Lc (sketched in Fig. |2J) . The re- 
sulting exact value for \t^(Ep, U)\ 2 deviates from the one 
obtained assuming the combination law of one-body scat- 
terers in series. This [/-dependent deviation is due to in- 
duced correlations in the attached leads, and its depen- 
dence on Lq allows to determine the size of the region 
which acts as an effective one-body scatterer. 



2 Embedding method, extrapolation and 
correlation cloud 

To study the mapping of a bare many-body scatterer cou- 
pled to leads onto an effective one-body scatterer with 
interaction-dependent coefficients, we take a model of N 
spinlcss fcrmions in a chain of L = Ls + Ll sites. The 
Hamiltonian (with even Ll) reads 



i=2 

Ll/2+Ls 
J=L L /2+2 



(1) 



The hopping amplitude tb. = 1 between nearest neigh- 
bor sites sets the energy scale, (cj) is the annihilation 

(creation) operator at site i, and n$ = cjcj. The near- 
est neighbor repulsion U acts upon Ls consecutive sites 
and gives rise to many-body scattering. We take a half- 
filled model (N = L/2), with a potential V+ — 1/2 being 
due to a positive background charge which exactly com- 
pensates the repulsion U inside the scatterer. Therefore, 
our model exhibits particle-hole symmetry and a uniform 
density, without Friedel oscillations around the scattering 
region where the fermions interact. 



The scattering geometry corresponds to two leads of 
Ll/2 — > oo sites connected by an interacting scatterer 
of Ls sites. The electrons do not interact in the leads, 
a necessary condition for having appropriate asymptotic 
scattering channels in one dimension. For the embedding 
method, we consider the ring geometry sketched in Fig.^ 
the scatterer being closed on itself via a non-interacting 
lead of Ll sites. This is achieved by adding a hopping term 

— thc\c L exp(i(f>) + h.c. 

to the Hamiltonian l|T)l. the flux 4> driving a persistent 
current J(U) in the ring. As the flux dependence of J(U) 
extrapolated to the limit Ll — > oo demonstrates UJ, the 
many-body scatterer behaves as an effective one-body 
scatterer, but with an interaction-dependent elastic trans- 
mission coefficient \t(Ep, U)\ 2 . Instead of using J(U), it is 
simpler [H] to get \t(Ep, U)\ 2 from the charge stiffness 

D(U, L S , L) = (-1)"- (E Q (U, L S , L) - E„(U, L s , L)) , 

(2) 

where Eq(U, Lg, L) — E„(U,Ls,L) is the change of the 
ground-state energy from periodic to antiperiodic bound- 
ary conditions. D(U, Ls, L) is obtained by the DMRG im- 
plementation for real Hamiltonians, which can be used to 
study with a great accuracy systems as large as L = 120 
sites with N = 60 particles. In the limit Ll — > oo, one 
gets the modulus 



\t(E F ,U)\ 



7T floo^Ls) 

2D oo (*7 = 0,L s ) 



(3) 



of the transmission amplitude through the scatterer of Ls 
sites, Doq(U — 0, Ls) being the charge stiffness of the 
same ring for [7 = 0. 

To take the limit Ll — » oo is one of the key points 
of the embedding method. This extrapolation is also re- 
quired for pure one-body scattering, where the finite size 
corrections to formula J3J can be expanded in powers 
of 1/L. For many-body scattering, the DMRG study gives 
an empirical scaling law 



D(U, L s ,L) = D 00 (U,L s )exp 



C(U,L S ) 



(4) 



obtained for large Ll and small Ls , which allows to deter- 
mine the asymptotic value D OQ (U 1 Ls) necessary to obtain 
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\t(E-p, U)\ . Expanding this scaling law gives 

D(U,L S ,L)- £UEW« B{U ' Ls) (5) 



when L is large enough, where 



B(U,Ls) = C(U,L s )D 00 (U,L a ) 



(6) 



This is a power law decay, and not an exponential de- 
cay with a characteristic scale above which the finite 
size correction can be neglected. Numerical data show 
that B(U, Lg) is important for intermediate interaction 
strengths U. But in the limits U — > (no scatter- 
ing, total transmission) and U — ► oo (total reflection, 
Doo{U — > oo, Ls) — > 0), the finite size corrections vanish 
and B(U,L S ) -> 0. 

However, since adding one-body potentials in the re- 
gion of the Ls sites yields a finite size correction to 
D(U, Ls, L) even when U — » 0, the interpretation of these 
corrections is not straightforward. They do not depend 
only on the correlations induced in the attached lead by 
the interaction acting inside the scatterer, but also on 
more trivial one-body aspects. 



3 Combination of two many-body scatterers 
in series 

A more direct approach, where the finite size effects are 
only due to many-body correlations, consists in taking two 
identical scatterers connected by a scattering-free lead of 
size Lc in which the electrons do not interact, as sketched 
in Fig. |21 Since the scattering channels begin at the first 
attached sites of the leads when U = 0, there is a sim- 
ple combination law for one-body scatterers in series. To 
study how this combination law is broken with increasing 
U when Lq is small allows to show that the size of the 
effective elastic scatterer is larger than the region where 
the carriers interact. When U ^ 0, the scattering chan- 
nels begin only asymptotically far from the many-body 
scatterer. 

Without interaction, a scatterer can be described at 
energy E-p by a unitary scattering matrix Ss, written in 
terms of its reflection and transmission amplitudes r, r' 
and t, t' as 

*-(;?)■ <7) 

The scatterer being symmetric upon time reversal, one 
has t — t' , while r = r' if the scatterer is symmetric upon 
space inversion. The transfer matrix Ms (giving the flux 
amplitudes at the right side in terms of the flux amplitudes 
at the left side) reads 



M s 



1/t r/t 
r* ft* 1/t* 



(8) 



The total Hamiltonian and the parity operator can be si- 
multaneously diagonalized if one has inversion symmetry, 
to give even and odd standing-wave solutions which can 
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Fig. 3. Conductance g(Lc) of the set-up sketched in Fig. 
with Ls = 2, for U = 1. The points are obtained with the 
embedding method for the ensemble. The dashed line at g — 
0.717 gives the approximate value obtained from Eq. 11511 . with 
\t\ obtained by the embedding method for a single scatterer. 
The solid line is a fit with the form g(L c ) = 0. 7174-0. 057/Lc • 
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Fig. 4. Conductance g(Lc) as in Fig. but for U = 2. The 
dashed line at g = 0.357 represents the value yielded by Eq. 
pSfr. and the solid line is the fit g(Lc) = 0.3572 - 0.077/Lc. 



be written as ip® — cos(ki + S ) and ipl = sin(ki + Si) at 
the right side of the scatterer, and t/>? = cos(fci — So) and 
ipl = s'm(ki — Si) at its left side. The two phase shifts 8q 
and Si are related |15| to t and r by 




t = (cxp(2i(5o) 
r = (cxp(2i(5o) 



exp(2W 1 ))/2, 
exp(2££i))/2 . 



(9) 



Due to symmetries, Ss or Ms have only two free pa- 
rameters: the modulus \t\ — cos(<5o — Si) and the phase 
a = 5q + Si of the transmission amplitude t, the unitar- 
ity of Ss (\t\ 2 + |r| 2 = 1 and r/r* = —t/t*) giving r. We 
can determine \t\ by the embedding method. The Friedel 
sum rule gives a. If one introduces a scatterer with 
inversion symmetry in the central region of a scattering 
free lead, this rule states f^Jj that 



a = 5q + Si = ttN{ 



(10) 



for spinless fermions in one dimension. Nf is the number 
of displaced fermions when the scatterer is introduced in 
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Fig. 5. Error Sg made when using the combination law of Eq. 
1141 for having the conductance of the set-up sketched in Fig. El 
with Ls = 2 and different values of U, as a function of an even 
number Lc of connecting sites. The lines give an A(U,Ls)/Lc 
fit. 



the central region. For a uniform filling factor v — 1/2, 
Nf = is/2, and the phase a reads 



a = irNf = = fcpLg . 



(11) 



where fcp = tt/2 is the Fermi wave number. For the spin- 
less case in one dimension with a uniform density, this 
simply means that the transmitted wave has Nf changes 
of sign when one transfers a fermion through a scatterer 
containing Nf others. This is obvious for U = as well 
as for (7^0. Using the same rule, the ideal ballistic 
lead of Lc sites has a modulus |t(Lc)| = 1 and a phase 
a(Lc) = kpLc- Its transfer matrix reads 



M fi = 









,ik F Lc 



(12) 



The combination law of one-body scatterers in series being 
a simple matrix multiplication for the transfer matrices, 
the total transfer matrix Mt(Ef) of the ensemble is given 

by 

M T = M s ■ M ■ M s . (13) 

For the total transmission coefficient | 2 through the en- 
semble, expressed in terms of the transmission t of each 
element and of Lc, this gives 



2(1 - |t| 2 )(l + cos (2fc F L c - 2a)) + \t\ 4 



(14) 



Since \t\ = 1 when is is odd [STITH] , we consider only 
even values of Ls- Taking a = ttLs/2 gives the Landauer 
conductance g = \t^\ 2 — 1 if Lc is odd and 



atp-2r 



(15) 




Fig. 6. Amplitude A(U, Ls) of the fits shown in Fig. as a 
function of U for different values of Ls- 



4 Correlation-induced deviations from the 
non-interacting combination law 

Fig. and ri] show the conductance g for two scatterers of 
Ls = 2 sites in series, as a function of the length Lc of 
the coupling lead. The data points are directly obtained 
from the embedding method, without assuming a combi- 
nation law for scatterers in series. Resonances with g = 1 
occur for odd Lc = {1, 3, 5, . . . }. For even Lc, the dashed 
lines represent the Lc-independent values |i| 4 /(|£| 2 — 2) 2 
implied by Eq. i|14[l . the coefficient \t\ being obtained us- 
ing the embedding method for a single scatterer. Within 
the accuracy of the extrapolation procedures required for 
having the transmission \t\ of an individual scatterer and 
the total conductance g, the result of (|14|l gives the cor- 
rect value when Lc — » oo, but overestimates g for small 
even values of Lc- The difference 



Sg{L c ) = g{L c 



9{Lc) 



(16) 



is shown in Fig.0for even Lc at different values of U. For 
even Lc, Sg(Lc) decays as a function of Lc as 



5g(L G ) 



A(U,L S ) 



(17) 



for Lc even. 



with an amplitude A(U, Ls) which is shown in Fig. as a 
function of the interaction strength U. This 1/Lc decay is 
reminiscent of the 1/L decay characterizing D(U 1 Ls, L) — 
DaoiU, Lg) and of the screening at large distances (larger 
than the Thomas- Fermi screening length) of the potential 
of a point charge by non-interacting electrons (Friedel os- 
cillations, RKKY interactions . . . ) in one dimension. This 
suggests that the decay could be faster for leads of higher 
dimensions (1/L C decay in d dimensions). The amplitude 
A(U, Ls) — > when U — * (one-body scatterers) and 
when U — > oo. In this latter limit, the scatterers become 
decoupled from the leads, the energy for an electron to 
enter or to leave a scatterer being oc U. A(U) is maxi- 
mum near U = 2, a value where in the thermodynamic 
limit Ls — > oo the Luttinger liquid becomes f^Hl a Mott 
insulator for spinless fermions. 
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For all odd values of Lq, the data for the total trans- 
mission coincide with the value g — 1 obtained from Eq. 
lfT4*|) . and Sg(Lc) = 0. The even-odd dependence on the 
parity of Lq shows that the convergences of the phase a 
and the modulus \t\ of the effective scatterer are character- 
ized by different scales. One has a — ttN{ across the scat- 
terer, directly on a scale L§, independently of U, while \t\ 
reaches its asymptotic value on a much larger scale. This is 
not surprising since a depends on the mean density, while 
|t| depends on the correlations of its fluctuations. In our 
model with a compensating background charge, the mean 
density does not exhibit Friedel oscillations. Let us under- 
line that the correlation clouds which have to be included 
with the many-body scatterer to form the effective one- 
body scatterer must not be confused with the screening 
clouds characterizing the charge density. 

5 Discussion of the relation to the Kondo 
problem 

To obtain the effective one-body matrix S(E-p) of a cor- 
related system of spinless fermions is a problem which 
displays a certain similarity with the Kondo problem of 
a spin degree of freedom surrounded by a metallic host. 
In the two cases, it is crucial to couple the many-body 
system to non-interacting conduction electrons. For the 
Kondo problem, the original 3d model can be mapped 
onto a Id lattice without interaction embedding a Hub- 
bard impurity. In Wilson's renormalization group trans- 
formations |12| . the embedded many-body Hamiltonian is 
progressively mapped onto an effective one-body Hamilto- 
nian describing the low energy states. In this transforma- 
tion, the coupling between different length scales is taken 
into account progressively, working out from the impu- 
rity to the longer length scales and lower energies. The 
states at sites near the impurity involve conduction states 
spanning the full band width 2£h, while the states located 
far from the impurity involve conduction states near the 
Fermi level, with a progressively reduced band width. This 
NRG method has been used recently [21] to calculate the 
effective one-body Hamiltonian of a few Hubbard sites em- 
bedded in a non-interacting chain, and the corresponding 
phase shifts. 

In our spinless case, the embedded system would give 
rise to inelastic and elastic scattering if it was in the vac- 
uum. Due to the attached non-interacting leads inelastic 
processes become progressively blocked by a Fermi vac- 
uum which eventually takes place in the leads, at a large 
distance from the scatterer. Our DMRG study leads us 
to a similar conclusion as NRG studies of impurities with 
spin. It shows that the effective one-body elastic scatterer 
necessary for extending the Landauer formulation of co- 
herent transport to correlated fermion systems must in- 
clude parts of the attached leads where the interacting re- 
gion induces power law correlations. Physically, this gives 
a dependence of the conductance of a mesoscopic scat- 
terer upon the nature of the used leads which depends on 
the strength of the interactions inside the scatterer. Even- 
tually, let us mention that the vertex corrections due to 



inelastic scattering vanish |21j when T — > 0, in a perturba- 
tive approach to the Kubo conductance of an interacting 
region embedded between semi-infinite leads. While this 
agrees with our findings for Ll — > oo, it is likely that these 
corrections do not vanish when Ll is finite, and exhibit 
similar power-law decays as Ll increases. 
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